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We present simple new approximate formulas, for both scalar and spinor QED, for the number of 
particles produced from vacuum by a time dependent electric field, incorporating the interference 
effects that arise from an arbitrary number of distinct semiclassical turning points. Such interference 
effects are important when the temporal profile of the laser pulse has subcycle structure. We show 
how the resulting semiclassical intuition may be used to guide the design of temporal profiles that 
enhance the momentum spectrum due to interference effects. The result is easy to implement and 
generally applicable to time-dependent tunneling problems, such as appear in many other contexts 
in particle and nuclear physics, condensed matter physics, atomic physics, chemical physics, and 
gravitational physics. 
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I. INTRODUCTION 

The Heisenberg-Schwinger effect is the non-perturbative production of electron-positron pairs when an external 
electric field is applied to the quantum electrodynamical (QED) vacuum [IHZ]- It was one of the first non-trivial 
predictions of QED, but the effect is so weak that it has not yet been directly observed. However, new experimental 
developments in ultra-high intensity lasers [SI (Sj may soon bring us to the verge of this extreme ultra-relativistic 
regime |10j . This experimental progress has renewed theoretical interest [llj . and recent results suggests that the 
effect may become observable in the 10^^ — lO^^W/cm^ intensity range, three or four orders of magnitude below the 
"Schwinger limit" of 4 x lO^^W/cm^ , which comes from an estimate based on a constant electric field. New theoretical 
ideas involve combining multiple copies of identical pulses [TH [13] , and also shaping pulses in special ways using the 
"dynamically assisted Schwinger mechanism" |14| . in which a superposition of two time-dependent pulses, one strong 
but slow, and the other weak but fast, can lead to a significant enhancement of the tunneling process associated with 
the Schwinger effect. An explicit experimental realization has been proposed [TS] that suggests an observable rate 
of particle production. A closely related theoretical idea is that of a "catalyzed Schwinger mechanism" [TB], which 
can also be viewed as photon-stimulated pair- production |T7] . realizing the more general mechanism of an induced 
metastable decay process [18]. The importance of cascading effects has been emphasized in [19]. These, and other 
theoretical analyses of more realistic laser fields, such as plane waves of finite extent [W, show that the precise form 
of the laser field can have a significant effect on the resulting pair production yield and momentum distribution. The 
strong sensitivity is not so surprising since it is a non-perturbative effect, but this makes it correspondingly difhcult 
to do precise computations. 

In the quantum field theoretic approach [2J, the theoretical problem is to compute the non-perturbative imaginary 
part of the "effective action", r[A] = h\n det [ip— m], where the Dirac operator, p= ^^{d^ — i^A^)^ defines the 
coupling between electrons and the applied (classical) electromagnetic field that represents the field produced by 
the laser pulse. However, computationally we are currently limited to one-dimensional fields such as time dependent 
electric fields -E(t), in which case the problem can be more conveniently expressed as a "quantum mechanical" 
scattering problem, invoking Feynman's picture of anti-particles as particles traveling backward in time 21 . This 
requires the computation of a reflection probability for an over-the-barrier scattering problem, that can be done 
numerically or using WKB [52H2S], or in the quantum kinetic approach [^fHSU] . The WKB approach is based on a 
relativistic extension of Keldysh's seminal work for atomic ionization in time-dependent electric fields 31]. Recently 
it has become clear that this WKB analysis must be extended to incorporate interference effects when the temporal 
profile of the laser pulse has sub-cycle structure. 

For example, it has been shown |32| that the momentum distribution of the produced particles is extremely sensitive 
to the "carrier phase" of a laser pulse, the phase offset between the pulse envelope and its oscillatory function. 
Moreover, this sensitivity reveals a distinct difference between spinor and scalar QED, which are conventionally treated 
on an equal footing at leading non-perturbative order. The oscillatory behavior of the momentum distribution is an 
interference effect, and can be quantitatively explained by incorporating the interference between different semiclassical 
saddle points. This was done in (33) for the case of two distinct saddle points, using the phase integral method and 
the Stokes phenomenon. In this paper. Here we use another even simpler method, based on the Riccati form of the 
scattering problem, and present new results for the case of an arbitrary number of distinct saddle points. 

Physically, such interference phenomena are familiar from strong-field atomic and molecular physics, discussed 
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long ago in the theory of atomic ionization , and observed experimentaUy in photoionization spectra [371 131] ■ 

These ideas have even led to the proposal for an all-optical double-slit experiment in the time domain, using vacuum 
polarization effects [511 132]- They also appear naturally in any time-dependent tunneling effect, such as the Landau- 
Zener effect or other condensed matter systems [ID] , chemical physics [H] , as well as gravitational H5] and particle 
physics [Hj. 

In Section II we briefly review the scattering formalism for the QED pair production effect. In Section III we 
recall the numerical approach, and present our approximate expressions for the particle number. Sections IV and V 
contain explicit examples of particular temporal profiles for the electric field E{t) that illustrate various features of 
the interference phenomena, and the final Section contains our conclusions. 



II. SCATTERING FORMALISM 



In this Section we recall briefly the scattering formalism of the pair production problem for both scalar and 
spinor QED, as we wish to compare the two cases in subsequent sections. For a linearly polarized electric field 
E = (0, 0, E{t)) that is time dependent and pointing in the direction, we choose a vector potential A = (0, 0, A{t)), 
with E{t) = —A{t). For such a field, spatial momentum is a good quantum number for the produced particles, so 
we can decompose the quantum field operators in terms of spatial momenta. For both scalar and spinor QED, the 
number of particles produced in each momentum mode can be expressed in terms of the reflection coefficient for an 
effective Schrodinger-like scattering problem. Physically, this is due to Feynman's interpretation of antiparticles as 
particles propagating backwards in time fSP , and has been used as a basic tool in the WKB analysis of the particle 
production problem [22.- 24^ i26j . The point of this current paper is to extend such semiclassical analyses to incorporate 
interference effects due to multiple saddle points, as this phenomenon naturally occurs for time dependent electric 
fields with sub-cycle structure, as is the case for more realistic representations of intense laser pulses. 



A. Scalar QED 



We decompose the scalar field operator as 

<I>(f , i) = / d'fc e'^-^ [MtW + Kit)blk) (1) 



where Ok and 61^. satisfy standard bosonic commutation relations, for each mode k. The Klein-Gordon equation for 
^{x^t) translates into the following equation for the mode functions 4>k{t): 

Mt) + QUt)Mt) = (2) 

where we define 

Ql{t) ^m^ + kl + (fell - qA{t)f (3) 
Equation ([2| has the form of a Schrodinger-like equation in the variable t 

- Mt) - - qA(t)fMt) - (™' + kl)Mt) (4) 

with "potential" V{t) = — (fc|| — qA{t))'^, and "energy" (m^ + k'^). We implement the Bogoliubov transformation by 
defining ak{t) and (3w{t) as follows: 

Mt) = -^Q^{t) ( -^d^e-^rQ. _ Mt^^.^Q^ 

\V^Q^) V^&) J 

For each mode, the Bogoliubov coefficients, at and /3k, satisfy the first-order coupled equations: 
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This Bogoliubov transformation implements a change from the time-independent basis of creation and annihilation 

-k' 

linear transformation: 



operators, Ok and ^l^., to a time-dependent basis of creation and annihilation operators, ak(i) and b^_j^{t), via the 



flk 
fo-k 



(8) 



The bosonic commutation relations are preserved by the unitarity condition: |Q!k(^)P ~ |/3k(i)P — 1- The number 
of pairs produced in the momentum mode k, from vacuum, is given in terms of the modulus of the coefficient /3k at 

t — -foo: 

Nk = \Mt = +cx^)\^ (9) 
The relation to quantum mechanical scattering arises because we can express -/Vk in terms of the reflection probability. 



l^kl^ 



t — -\-OQ 



, for the effective " Schrodinger" problem ( 2 1 



l^kl 



Recall from Q that this describes the situation of over-the-barrier scattering, so the reflection probability is expo- 
nentially small, and so we can often make the approximation: iVk ~ |iikP- 



B. Spinor QED 

An analogous mode decomposition exists for spinor QED. We expand the spinor field operator '^{x,t) as: 

^{x, t) = Y. f ^'" ^ "k.s + l^-k,s(i) folk.s) (11) 

where Ok and 6^^. satisfy standard fermionic anti-commutation relations, for each mode k, and the sum is over helicity 
s = ±1. In a suitable Dirac matrix basis, the time dependent spinors, Uk,s(i) and Wk,s(i)j can be written in terms of 
a single complex function ■(/'k(i) that satisfies the Schrodinger-like equation: 

V^k (i) + (Qk W + ^-^ii it)) Mt) = (12) 

We implement the Bogoliubov transformation by defining ai^{t) and /?k(0 follows: 

^..(t)= e-'^'q--f ^"^^^^ ^'^'g- (13) 

V2Qk(i)(Qk(t)-fc||(i)) V2Qk(t)(Qk(t) + A:||(t)) 

^u(t) = -tOu(t) ( "'-(^) .-»rQ. _ ^k(t) ,^rQ.\ 

' ' l^V2Qk(t)(Qk(i)-A:j|(i)) v/2Qk(i)(Qk(i) + fci!(i)) J 

For each mode, the Bogoliubov coefficients satisfy the first-order coupled equations: 

M) = -i&ak(()e-"-r'«'. (15) 



where = m? + k'j_. This Bogoliubov transformation implements a change from the time-independent basis of 
creation and annihilation operators to a time-dependent basis of creation and annihilation operators, with the unitarity 
condition: |Q!k(i)P + |/3k(^)P = 1- Note the opposite sign from the scalar QED case. The number of pairs produced 
in the momentum mode k, from vacuum, is given in terms of the modulus of the coefficient /3k at t = +oo: 

Wk = |/3k(i = +oo)|2 



where |i?kP is defined as the reflection probability, |-RkP 
so we often make the approximation: A'k ~ |Ji!kP- 



l + |i?kP 

2 



(16) 



/3k (t) 



. Again, the reflection is typically very small. 
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III. SCATTERING FORMALISM: NUMERICAL AND SEMICLASSICAL APPROACHES 

A. Scalar QED: numerical computation 

It is straightforward to convert the Schrodinger-Uke scattering problem Q into a Riccati equation [53] that is 
suitable for simple numerical evaluation. (Wc now suppress the momentum mode label, k, since all modes decouple, 
and so can be treated separately.) From equations Q it is clear that the reflection amplitude, R — /3/a, evolves with 
time as: 

• a/S — /3 a 
li = ^ 

This Riccati equation is trivial to integrate numerically, for a given A{t) and longitudinal momentum kn , with the 



initial condition R{—oo) = 0, to obtain R{oo), whose magnitude squared gives the particle number (10). As discussed 
in |45| . this is completely equivalent to the quantum kinetic equation approach. We will use this numerical formalism 
in order to obtain "exact" particle spectra, with which we can compare our semiclassical approximations. 

B. Spinor QED: numerical computation 



For spinor QED, the argument is very similar. From the equations (15) it is clear that the reflection amplitude, 
R — f3/a, evolves with time as: 



^ ^ a$ - fig 



fc|| e± 



Note the different signs from the scalar case (17), and the different form of the function out the front. Again, it is 



simple to implement numerically, with the initial condition R(—oo) = 0, to obtain i?(+oo), and hence the particle 



number from (16). 

C. Scalar QED: semiclassical approximation 



To motivate a semiclassical approximation to the Riccati equation (17), consider the fact that R{t) is always small. 



and so neglect the nonlinear term on the right-hand-side |46j . Then we have simply, 

i?(oo) ~ P ^ /-oo Q(*') d^' dt (19) 

This integral is dominated by the contributions of the poles, where Q = 0, which are the semiclassical turning points 
tp. In the neighborhood of such a turning point, change variables from t to the "singulant" function 

m = J Q{t') dt' (20) 

Now assume there is a first order zero of (t) [a similar argument applies for other orders of poles] , so that near the 
turning point, Q ^ c^t — tp, and ^ ~ ^c{t~ tp)^^"^ + ^p. Then the approximate equation for R(t) can be expressed 
as 

(21) 

Therefore, each pole will contribute a term i?(oo) ss — ^e"^*^'', where = fy^Q{t) dt, and we have chosen 
to refer all the phase integrals to t = — cx). In fact, this approximation does not give the correct pref actor. This 
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prefactor problem was noted already in the seminal papers [17], where it was resolved by comparison with soluble 



cases. In order to obtain the correct prefactor, we must consider also the nonlinear term in (17), and keep all the 
multiple- integral iteration terms. This procedure yields a prefactor of magnitude 1 Using these results, we 

obtain an approximate expression for the reflection amplitude with a contribution from each turning point in the 
upper half complex plane; 

i?(oo)«^e-2^^-"-'3W<it (22) 

Integrals of Q(t) along the real axis are real, while those along the imaginary direction are imaginary, so it is natural 
to split the exponents into phases and real parts. Let us define Sp = Re(ip) as the real part of a complex turning 

point tp. Then we can separate out a common phase factor e~'^^-^-'=° in the sum in (22), and write 



where the phase, 9p — J^^ Q{t) dt, is the phase accumulated by integrating Q{t) along the real axis between neighboring 
turning points. These phases incorporate the interference effect between distinct turning points and yield a simple 



expression for the reflection probability, when we take the modulus squared of the reflection amplitude in ( 23 ) 



Within this approximation, there is actually no distinction between the particle number N in ( 10 ) and the reflection 
probability so we obtain the approximate expression: 

^.calar ^ ^ g-2<'" ^ ^ 2 COS (2 f?'^'"'') g-^^''-^"'' (24) 



where we have defined 



4'^ ^ I / Q^{t)dt\ (25) 



,{t)dt (26) 



We have restored the momentum label k to emphasize the fact that the answer depends on k, because (5k(0 depends 
on k. It should of course be remembered that this means that the location of the turning points tp also depends on 



k, and so do the interference terms 9^^'^ \ The first term in (|24|) is the sum over the contributions of independent 



turning points, while the second sum characterizes the interference between different turning points. The dominant 

(p) 

contributions are from turning points with the smallest values of , and interference effects are significant for pairs 
of turning points for which these integrals are comparable in magnitude. Loosely speaking, this often corresponds 
to a rule of thumb that turning points closest to the real axis tend to dominate, and interference effcts are strongest 
between pairs of turning points that have comparable distance from the real axis. 

For later use, we record the approximate expressions for one, two and three complex conjugate pairs of turning 
points. If a single turning point pair dominates, then we have the familiar textbook expression |52) : 

^scalar _ ^-2kI^^ (27) 

If there are two pairs of turning points, (^1,^^^) and {t2,t2), with comparable real exponential factors e"^^^ and 
e~^^k ' ^ then there is a single interference term 

^scalar ^ g-2<^) _^ g-2if<^' ^ 3 COS (2 (28) 

where 6*^^'^^ = J'^ Qk{t) dt. This is the case that was studied in [33] ■ If there are three turning point pairs, {ti,tl), 
(t2,t2)i s-iid (isjig), each with comparable real exponential factors e~'^^'^\ then there are three interference terms: 

^scalar ^ g-2<> ^ ^-2<^) ^ ^-2<^) ^ 3 COS (29i'^'^) g-^^'-^^' 

+2 COS (29^^''^) e-<^'-<^' + 2cos (29i' '^) e~<''~<'' (29) 
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where 6'^^'^'' — J^^ Q'k(t)dt, o'^'^^ = J^^ Qkit) dt, and 0^'^^ = J^^ Q]s.{t)dt. The extension to more pairs of turning 
points is clear. 

In the next Section we will illustrate these interference effects with explicit examples of electric fields that produce 
exactly one, two and three pairs of turning points. In the semiclassical regime, the expression (24 1 is an excellent 
approximation, and describes the interference effects both qualitatively and quantitatively for a broad range of physical 
parameters. 



D. Spinor QED: semiclassical approximation 



For spinor QED, we can apply a similar argument to the Riccati equation (18). The difference is that the initial 
approximation, which is then iterated, yields a different function: 

Ji^^'^'^^dt (30) 

In the vicinity of a turning point, we have 

2Q2 2Q fc - A{t) 2Q ^ ' 

with the sign depending on the branch. These signs alternate between successive turning points, so we obtain an 
extra (alternating sign) phase 

R{oo) « ^(-1)^' e''^/^ Qit)dt (32) 

which leads to an approximate expression for the particle number for spinor QED: 

^,-2<-) + ^ 2cos (2ei^'P">) {-1)(P-P') (33) 



^spinor 



where K^^ and 9^'^ -* are defined exactly as in (25 1 and (26). The only difference from the scalar QED case lies in 
the signs of the interference terms. For example, if a single turning point pair dominates, because of a dominant real 

factor e-^' ^1 , then there is no interference and we have just 

^^Pinor _ (34) 



which is the same as for scalar QED. If there are two pairs of turning points, (ti,tl) and (^2,^2)1 with comparable 

real exponei 
scalar case: 



real exponential factors e ' and e ' ^ then there is a single interference term, with the opposite sign from the 



This is the case that was studied in i33j- If there are three turning point pairs, (ii,i|), (i2j^2)i ^^"^ (^Sj^s); each with 
comparable real exponential factors e^'^^i^' , then there are three interference terms, with signs as follows: 

^ e-2<%e-^<%e-^<^'-2cos(2 0(^'^))e-<'-<^' 

-2cos {29^^''^) e-<^'-<^' + 2C0S {29^'^) e-<''~<'' (36) 

where 9^^'^^ = J^^ Q\^{t)dt, 9^^'^^ — J^^ Qkit) dt, and 9^'^"^ — J^^ Q\^{t)dt. The extension to more pairs of turning 
points is clear. 



IV. ILLUSTRATIVE EXAMPLES 



In this section, we compare our semiclassical approximations (24 1 and (33) that incorporate interference effects, 
with the (exact) numerical approach based on the Riccati equations ( 17 ) and ( 18 ) , for both scalar and spinor QED. For 
this comparison, we have constructed electric fields such that the corresponding over-the-barrier scattering problem 
has precisely one, two and three pairs of complex conjugate turning points. 
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A. One Pair of Turning Points 



An example of a gauge field with only single pair of turning points is the single-bump field 

Eq 



m = 



(37) 



where Eq is the field strength amplitude, and ui is the inverse width, as shown in the left panel of Figure [T] The 
associated vector potential can be taken as 



Ait) 



Eot 



(38) 



This vector potential is plotted in the right panel of Figure [l] Note that E{t) is an even function, while A{t) is an 




t 




FIG. 1: The form of the electric field E{t) in ( 37 1, and corresponding vector potential A{t) in (38 1 for a single complex conjugate 



pair of turning points. E{t) is an even function, while A{t) is an odd function. 

odd function. This field has exactly one pair of complex conjugate turning points, (ti(fc), tj(fc)), with 

-(fc-z) 



hik) 



^El + a;2 + 2ikufl - k^uj^ 



(39) 



Note that as a function of the longitudinal momentum, k, the pair of turning points moves around in the complex 
plane, as shown in Figure [2j but remain a complex conjugate pair. 




FIG. 2: The locations of the complex conjugate pair of turning points, in the complex t plane, for three different values of 
longitudinal momentum. These plots are for the vector potential A(t) in 1 38 1, with Eq — 0.1 and cj = 0.1, for longitudinal 
momentum vales fcy — —1 (left), fcy = (center), and fcy = 1 (right), in units with m = 1. Note that the turning points are 
closest to the real axis for k\\ = 0. 



Figure 3 shows a comparison between the approximations (27) and (34) and the exact numerical results, for the 
particle number as a function of longitudinal momentum. The left plot is for scalar QED and the right plot is for 
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spinor QED. There is no oscillatory structure in this momentum spectrum, as expected since there is no interference 
term for just a single pair of turning points. Thus, the spectra are the same for scalar and spinor QED, and there 
is good agreement between the approximate and exact results. While the form of the electric field (37) was chosen 
so that there is precisely one complex conjugate pair of turning points, similar behavior is obtained for other "single- 
bump" electric fields such as E{t) — Eq sech^ (ujt) , or E{t) — Eoexp^—uj^t^), for which there is an infinite tower of 
turning points pairs, but only one pair [the one closest to the real axis] dominates, and the approximate expressions 
(27) and (34) again provide extremely accurate answers. 




FIG. 3: Scalar (left) and spinor (right) QED momentum spectra for vacuum pair production, as a function of longitudinal 
momentum, for the electric field (37 1 that has one pair of turni ng p oints. Th e thick (blue) lines show the numerical calculation, 
and the dashed (red) lines show the approximate expressions (27 1 and (34 1. The field parameters were chosen as: Eq = 0.1, 
and ij — 0.1, in units with m — 1. 



B. Two Pairs of Turning Points 



To illustrate the effect of interference between pairs of turning points, we now consider an example of a vector 
potential leading to precisely two pairs of complex conjugate turning points. This field was considered already in |33j . 
and here we give more details. Consider the electric field 

Eit) = 2^°^ (40) 

where Eq is the field strength amplitude, and uj is the inverse width, as shown in the left panel of Figure |4j The 
associated vector potential can be taken as 



This vector potential is plotted in the right panel of Figure [4] Note that E{t) is an odd function, while A{t) is an 




FIG. 4: The form of the electric field E{t) in (40 1, and corresponding vector potential A{t) in (41 1 for two complex conjugate 
pairs of turning points. E{t) is an odd function, while A{t) is an even function. 



9 



even function. For this vector potential, there are two complex conjugate pairs of turning points, (ti(fc), tj(fc)), and 
(t2(fc),t5(fc)), where: 



ti(fc) 
t2{k) 



\/—Eq — k(jj + iu! 



^3/2 

V^Eq - kuj - 



+ 1 



(42) 
(43) 



These turning points are illustrated in Figure [Sj An important difference from the case of a single pair of turning 
points shown in Figure [2] is that now the point of closest approach of the turning points to the real axis occurs at a 
nonzero value of fcy. This is reflected in the momentum spectrum for the two-pair case, shown in Figure 6, which is 
centered around a non-zero value of fc|| , while the momentum spectrum for the single-pair case, shown in Figure 3, is 
centered around /cy = 0. Also, observe that since the two pairs are equidistant from the real axis, we should expect 
strong interference effects between the two pairs of turning points, as indeed is seen in Figure |6] for both scalar and 
spinor QED. 



FIG. 5: The locations of the complex conjugate pair of turning points, in the complex t plane, for three different values of 
longitudinal momentum. These plots are for the vector potential A(t) in | |41[ ), with Eg — 0.1 and uj = 0.1, for longitudinal 
momentum vales fcy = (left), fcy — 1 (center), and fcy = 2 (right), in units with m — 1. Note that the two turning points are 
always equidistant from the real axis, and note that they are closest to the real axis at a nonzero value of fcy , which for these 
parameters is fc|| ~ 1.2. 



Figure |6] shows a comparison between the approximations (28) and (35) and the exact numerical results, for the 
particle number as a function of longitudinal momentum. Note the oscillatory behavior of the spectrum, due to 
the interference terms. Also notice that the interference term has the opposite sign for scalar and spinor QED, as 
reflected in the exact momentum spectrum. The agreement between the exact numerical results [solid, blue lines] and 
the approximate semiclassical expressions [dashed, red lines] is extremely good, both qualitatively and quantitatively. 

The form of the electric field ( 40 ) was chosen so that there are precisely two complex conjugate pairs of turning 
points. In fact, for other electric fields with temporal profile that is an odd function of t, as in Figure |4l we find 
that there are two dominant pairs of turning points. For example, this occurs when E{t) = EQUjtsech {ujt), or 
E{t) = EqujI exp(— i^), for which there is an infinite tower of turning points pairs, but only two pairs [those closest 
to the real axis] dominate, and the approximate expressions (28) and (35) again provide extremely accurate answers. 



C. Three Pairs of Turning Points 



To illustrate further the effect of interference between pairs of turning points, we now consider an example of a 
vector potential leading to precisely three pairs of complex conjugate turning points. This goes beyond the field 
considered already in |33j . and permits us to verify the sign pattern of the interference terms in the spinor QED 
expression (36). Consider the electric field 



E{t) = 



Eq (1 - (3^^ 



2u^ 

5/2 



(44) 
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FIG. 6: Scalar (left) and spinor (right) QED momentum spectra for vacuum pair production, as a function of longitudinal 
momentum, for the electric field (40 1 that has two pairs of turning points. The thick (blue) lines show the numerical calculation, 
and the dashed (red) lines show the approximate expressions ( |28[ ) and ( |35[ ). The field parameters were chosen as: -Eo = 0.1, 
and iu — 0.1, in units with m — 1. 



where Eq is the field strength amplitude, and wi and uj2 represent two independent inverse width scales. The form of 
this electric field is shown in the left panel of Figure [7] The associated vector potential can be taken as 



Ait) 



Eot (l - 
" (l+a;2t2)3/2 



(45) 



This vector potential is plotted in the right panel of Figure [7j Note that E{t) is an even function, while A{t) is an 





FIG. 7: The form of the electric field E{t) in (44 1, and corresponding vector potential A{t) in (45 1 for three complex conjugate 
pairs of turning points. E{t) is an even function, while A{t) is an odd function. 



odd function, as in the case of one pair of turning points, shown in Figure [TJ The equation for the turning points 
is a cubic equation in t^, so we obtain three complex conjugate pairs of turning points, (ii(fc), t^(/c)), (<2(fc), i2(^))i 
and {t3{k),t'^{k)). We do not write the expressions explicitly, as they are long and not particularly instructive. These 
turning points are illustrated in Figure [8) for various values of the longitudinal momentum. As in the case of a single 
pair of turning points shown in Figure|2l the point of closest approach of the turning points to the real axis occurs 
at /cji = 0. This is refiected in the momentum spectrum for the three-pair case, shown in Figure 9, which is centered 
around fc|| = 0, but in contrast to the momentum spectrum for the two-pair case, shown in Figure 3, which is centered 
around a nonzero value of fcy . Also, observe that since the three pairs arc almost equidistant from the real axis, we 
should expect significant interference effects between the various pairs of turning points, as indeed is seen in Figure [9] 
for both scalar and spinor QED. The momentum spectrum resulting from vector potentials that are odd functions of 
time [and hence electric fields that are even functions of time] exhibit symmetric oscillations centered around fcu = 0, 
since for odd gauge fields the phase integrands have the symmetry: Qk^^ (t) = Q-k^ (—t)- Therefore equations ([7| and 
( 15 1 remain invariant under the transformations fc|| — >■ — fc|| , and t — > —t. In the WKB framework, this fact is manifest 
as the symmetry of the turning point distribution under fcy — fcy, as can be seen from Figure |8] 



Figure |9] shows a comparison between the approximations (29) and (36) and the exact numerical results, for the 
particle number as a function of longitudinal momentum. Note the oscillatory behavior of the spectrum, due to the 
interference terms. Also notice that the interference terms have different signs for scalar and spinor QED, leading to 
different oscillatory behavior in the longitudinal momentum spectrum. The agreement between the exact numerical 
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FIG. 8: The locations of the complex conjugate pair of turning points, i n th e complex t plane, for three different values of 
longitudinal momentum. These plots are for the vector potential A{t) in (45 I, with Eq = 0.1, oji = 0.1, and ijJ2 — 1/15, for 
longitudinal momentum vales fcy = —1 (left), fcy = (center), and fcy = 1 (right), 
suggests we should expect significant interference effects, especially near fcii — 0. 



in units with m = 1. Their distribution 




FIG. 9: Scalar (left) and spinor (right) QED momentum spectra for vacuum pair production, as a function of longitudinal 
momentum, for the electric field ( |44[ ) that has three pairs of turning points. The thick (blue) lines show the numerical calculation, 
and the dashed (red) lines show the approximate expressions (29 1 and (36 1. The field parameters were chosen as: Eq = 0.1, 
oji — 0.1, and uj2 = 1/15, in units with m = 1. 



V. PULSE CONFIGURATIONS WITH FLAT ENVELOPES 



A significant advantage of the semiclassical approach is that it provides us with some physical intuition to guide the 
problem of designing the temporal shape of the electric field E{t) in order to produce a desired momentum specturm. 
This is an interesting, and difficult, "inverse problem", and in this Section we illustrate the idea with some examples. 
The treatment of temporally localized electric fields with sub-cycle structure is important within the context of 
vacuum pair production in that these type of fields represent more realistic pulse configurations with rich structure of 
momentum spectrum for the produced pairs. Further, from the experimental point of view, the investigation of field 
parameter dependance of the spectrum might be useful for achieving more prolific pair production. 

The first observation is that interference effects are more likely with an electric field with more temporal structure: 
the single-bump field exhibits no interference, while the fields with increasing numbers of maxima and minima tend to 
increase the level of interference in the momentum spectrum. This is obvious from the scattering picture, but it is not 
the whole story. We also see that there is a marked difference between cases where E{t) is an even or odd function. 
Indeed, for a complicated form [with many oscillations and possibly an envelope] of E(t), and hence correspondingly 
for A{t), the effective scattering potential — (fcy — A{t))^ changes dramatically as a function of fc|| , and it is not easy to 
see from the form of this scattering potential when there would be a minimum or maximum of the particle number. 
The best indicator comes from looking at the location of the turning points in the complex plane. This also shows us 
that interference effects will be most pronounced when different sets of turning points are approximately equidistant 
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from the real axis. 

We can illustrate these trends with some electric field configurations looking more and more like realistic laser 
pulses, with sub-cycle structure. In [35], the effect of the carrier phase was investigated for fields of a given frequency 
w, convolved with a Gaussian envelope function, with a phase offset (p: 

E{t) = Ea cos(w t + (/)) e-*'/(2^') (46) 

Strong interference effects are seen for the odd field where — 7r/2, and we now understand this as due to the 
interference between two dominant pairs of turning points |33j . Now we ask what happens if we change the shape 
of the field so that more than two pairs of turning points contribute. There should then be stronger interference 

4/4: S / S 

effects. This can be achieved by "flattening" the envelope function from a Gaussian to a factor e^* '"^ or e^* '"^ . 
We show below that this simple change in the envelope function increases the number of relevant turning point pairs, 
and correspondingly has a significant effect on the interference terms. We consider such envelope functions both for 
"cosine-like" and "sine-like" electric fields, corresponding to carrier phases (j) — and — 7r/2, respectively. 



A. Envelope Functions: exp[— f''/r*] 

Consider electric fields with an envelope function e~ ''^ , which is "flatter" than a Gaussian envelope. This leads 
to more turning points with approximately equal real parts, and therefore to stronger interference effects. Specifically, 
we first take an electric field temporal profile that is an even function of time 



Eoe- 



'^Lu cos(tcj) — At^ sin(tw)) 



(47) 



which comes from an odd vector potential 



A{t) = -Eo/cue-*"/'''' sm{ujt) 



(48) 



The forms of these fields are plotted in Figure [TOj and the turning point distribution is sketched in Figure 11 We see 
that there are more turning point pairs approximately equidistant from the real axis, suggesting stronger interference 
effects. The results for the produced particle number, as a function of longitudinal momentum, are shown in Figure 
12 for both scalar QED [solid, blue curve] and spinor QED [dashed, red curve]. Notice the single-peak structure for 
scalar QED, and double-peak structure for spinor QED, a reflection of the opposite sign of interference terms. 





FIG. 10; The form of the electric field E{t) in (471, and corresponding vector potential A{t) in (481 for two complex conjugate 
pairs of turning points. E{t) is an even function, while A{t) is an odd function. 



As a second example with the same envelope function, we consider an electric field temporal profile that is an odd 
function of time 



Eodd{t) 

which comes from a vector potential 



j4 

Eoe^^ (4i^cos(tw) + t^uj sm{tuj)) 



(49) 



A{t) = -Eo/cuer*^/^" cos(wi) 



(50) 
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FIG. 11: The locations of the complex conjugate pair of turning points, i n th e complex t plane, for three different values of 
longitudinal momentum. These plots are for the vector potential A{t) in (481, with Eo = 0.1 and uj = 0.1, for longitudinal 
momentum vales fci 



'1 (left), fell = (center), and fcy = 1 (right), in units with m = I. 



N 



7.x 10"' 
6.x 10"' 
5.x 10"' 
4.x 10"' 
3.x 10"' 
2.x 10"' 
1.x 10"' 




A(f)= -c"''/''sin((j) 



FIG. 12: The particle numbers for vacuum pair production, as a function of longitudinal momentum, for the electric field (47 1 



with the solid (blue) line showing scalar QED and the dashed (red) line showing spinor QED. The field parameters Eo, ui, and 
T were chosen as: Eo = 0.1, uj = 0.5, and r — 0.05, in units with m — 1. 



The forms of these fields are plotted in Figure [Tsj and the turning point distribution is sketched in Figure 14 We see 
that there are more turning point pairs approximately equidistant from the real axis, suggesting stronger interference 
effects. The results for the produced particle number, as a function of longitudinal momentum, are shown in Figure [15] 
for both scalar QED [solid, blue curve] and spinor QED [dashed, red curve]. Again, notice the single-peak structure 
for scalar QED, and double-peak structure for spinor QED, but now observe the asymmetry of the spinor spectrum. 





FIG. 13: The form of the electric field E(t) in ( |49[ ), and corresponding vector potential A{t) in ( |50[ ) for two complex conjugate 
pairs of turning points. E{t) is an odd function, while A{t) is an even function. 
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FIG. 14: The locations of the complex conjugate pair of turning points, i n th e complex t plane, for three different values of 
longitudinal momentum. These plots are for the vector potential A{t) in 1 50 1, with Eq — 0.1 and uj = 0.1, for longitudinal 



momentum vales fcy = — 1 (left), A;|| = (center), and fcy = 1 (right), in units with m = I. 



2.x 10-' 




A(t)= -e"'V^Cos(£j) 



FIG. 15: The particle numbers for vacuum pair production, as a function of longitudinal momentum, for the electric field (49 1, 
with the solid (blue) line showing scalar QED and the dashed (red) line showing spinor QED. The field parameters Eq, lo, and 
r were chosen as: Eo = 0.1, uj = 0.5, and r = 0.05, in units with m = 1. 



B. Envelope Functions: exp[— t**/r*] 

Now consider fields with an even flatter envelope function: e"***/"^ . This leads to even more turning points with 
approximately equal real parts, and therefore to even stronger interference effects. Specifically, we first take an electric 
field temporal profile that is an even function of time 

Ene ^ (t^uj cos(tu>) - 8f sm(tLu)) 
i^cvc„(t) = — ^ (51) 

which comes from an odd vector potential 

A{t) = -Ea/iue-*''/^^ sm{ujt) (52) 

The forms of these fields are plotted in Figure [l6j and the turning point distribution is sketched in Figure [iTj We 
see that there are even more turning point pairs approximately equidistant from the real axis, suggesting stronger 
interference effects. The results for the produced particle number, as a function of longitudinal momentum, are shown 
in Figure [Tsj for both scalar QED [solid, blue curve] and spinor QED [dashed, red curve]. Notice the very different 
forms of the momentum spectra, and in particular notice that the peak values for spinor QED are almost an order of 
magnitude greater than for scalar QED. 

As a second example with the same envelope function, we consider an electric field temporal profile that is an odd 
function of time 

-C'oddl.Ij = S [o6) 
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FIG. 16: The form of the electric field E{t) in (511, and corresponding vector potential A{t) in (52 1 for two complex conjugate 



pairs of turning points. E{t) is an even function, while A{t) is an odd function. 



• •••••• 



FIG. 17: The locations of the complex conjugate pair of turning points, in the complex t plane, for three different values of 
longitudinal momentum. These plots are for the vector potential A{t) in (52 1, with Eo = 0.1 and oj = 0.1, for longitudinal 



momentum vales fcy = — 1 (left), fc|| = (center), and — 1 (right), in units with m — 1 



which comes from a vector potential 



A{t) = -Eo/uje 



cos(aji) 



(54) 



The forms of these fields are plotted in Figure [19] and the turning point distribution is sketched in Figure [20] Again, 
we see that there are even more turning point pairs approximately equidistant from the real axis, suggesting stronger 
interference effects. The results for the produced particle number, as a function of longitudinal momentum, are shown 
in Figure 21 for both scalar QED [solid, blue curve] and spinor QED [dashed, red curve]. Notice the very different 
form of the spectra, and note that again the spinor QED peaks are noticeably higher than those for scalar QED. 



1. x 10"' 
8.x 10"" 
6.x 10"* 
4.x 10"* 

2. x 10"* 



A(t)= Sin ((J ) 



FIG. 18: The particle numbers for vacuum pair production, as a function of longitudinal momentum, for the electric field ( |51[ l, 
with the solid (blue) line showing scalar QED and the dashed (red) line showing spinor QED. The field parameters Eo, ui, and 
T were chosen as: Eq = 0.1, uj = 0.5, and r — 0.05, in units with m = 1. 
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FIG. 19: The form of the electric field E{t) in (531, and corresponding vector potential A{t) in (54 1 for two complex conjugate 



pairs of turning points. E{t) is an odd function, while A(t) is an even function. 



• • • • 



FIG. 20: The locations of the complex conjugate pair of turning points, in the complex t plane, for three different values of 
longitudinal momentum. These plots are for the vector potential A{t) in (54 1, with Eo = 0.1 and oj = 0.1, for longitudinal 



momentum vales fcy = — 1 (left), fc|| = (center), and fcy = 1 (right), in units with m = I. 



These examples clearly show that the flatter the envelope function, the stronger the interference effects, and with 
such a large number of turning points participating, there can be large differences between the pair production for 



spinor and scalar QED. For example, in Figure 18 we see almost an order of magnitude difference between the spinor 
and scalar QED peak particle numbers. Contrast this with the earlier examples, where even though the interference 
effects have different signs, they do not conspire to increase the overall magnitude of the peak values. 



N 



7.x 10-' 
6.x 10"* 
5.x 10-* 
4.x 10-" 
3.x 10" 
2.x 10" 
1.x 10" 




A(t)= -e-''A"cos(6i) 



FIG. 21: The particle numbers for vacuum pair production, as a function of longitudinal momentum, for the electric field (53 1, 
with the solid (blue) line showing scalar QED and the dashed (red) line showing spinor QED. The field parameters Eq, oj, and 
T were chosen as: Eq = 0.1, uj = 0.5, and r — 0.05, in units with m — 1. 
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VI. CONCLUSIONS 



In this paper we have investigated interference effects in the longitudinal momentum spectrum for particles produced 
from vacuum by a linearly polarized electric field that is spatially uniform but time dependent. The interference is 
due to the interaction between multiple semiclassical turning points, and becomes important when the temporal 
profiles have subcycle structure, as is true for more realistic laser pulse fields than jus t th e well -stu died single-bump 



fields hke E{t) = Eq sech^ {uj t) . We have given a simple new approximate formulas, fl24h and pSfc, for the number 
of produced particles, as a function of longitudinal momentum, for both scalar and spinor QED, for an arbitrary 
number of turning points, extending the result of |33| for the interference between two distinct turning points. As 
expected, the interference terms have different signs depending on the particle statistics. We have confirmed that 
these approximate expressions agree very well with the exact results, obtained by numerical integration of the Riccati 
form of the corresponding scattering problem, for electric fields having precisely one, two, and three complex conjugate 
pairs of semiclassical turning points. The approximate expressions provide important physical intuition that may be 
used to guide the shaping of the temporal profile of electric field pulses in order to obtain a particular momentum 
spectrum. In particular, we have shown that flattening the temporal envelope function leads to stronger interference 
effects, since more turning points interfere, and tends to increase the particle number for spinor QED relative to scalar 
QED. We hope that this semiclassical approach may be useful in guiding the design of planned laser experiments 
in order to observe this elusive non-perturbative Heisenberg-Schwinger effect for the first time. For example, the 
recent numerical results of Orthaber et al |53j concerning the momentum spectrum of vacuum particle production 
for the dynamically assisted Schwinger mechanism [TJj, in which a strong enhancement is seen when a weak but 
rapidly varying field is superimposed on a stronger but slower field, can be understood semiclassically in terms of 
the appearance of new saddle points that arise due to the additional weak field. In addition, such time-dependent 
tunneling problems appear in many other contexts [54], in particle and nuclear physics, condensed matter physics, 
atomic physics, chemical physics, and gravitational physics, and we anticipate that the simplicity of these results may 
prove useful in these other areas also. Finally, the semiclassical perspective in terms of interfering saddle points may 
prove useful in the search for a computationally effective formalism that also incorporates spatial inhomogeneities of 
the laser pulses, for example using worldline instantons ^55j or Wigner function methods |561 157j . 

We acknowledge support from the DOE grant DE-FG02-92ER40716. 
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